ON THE SPECTRAL GAP OF THE KAC WALK AND OTHER 
BINARY COLLISION PROCESSES 

PIETRO CAPUTO 

Abstract. We give a new and elementary computation of the spectral gap of the 
Kac walk on S . The result is obtained as a by-product of a more general obser- 
vation which allows to reduce the analysis of the spectral gap of an iV-component 
OQ ' system to that of the same system for N = 3. The method applies to a number of 

random "binary collision" processes with complete-graph structure, including non- 
homogeneous examples such as exclusion and colored exclusion processes with site 
disorder. 
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^— > . 1. Introduction 

The following model for energy preserving binary collisions was introduced by M. 
Kac in [8j. Let v denote the uniform probability measure on the sphere 

N 

> 5^ = 1^^: X> 2 = 1}, 

i ' »=i 



£ 



and consider the z^-reversible Markov process on S 1 with infinitesimal generator 
given by 

O " 1 N 1 r 2n r i 

cm = ^ E r / [fWrt ~ f{v) \ de ' (L1) 

where R l g , i ^ j is a clockwise rotation of angle 6 in the plane (r]i,rij)- As a convention, 
/\ we take R % g = Id. In words, the associated Markov process goes as follows: we have 

independent Poisson clocks of rate 1/2 at each coordinate; when coordinate i rings we 
choose uniformly at random (with replacement) another coordinate j; if j ^ i then we 
perform a rotation of angle 9 in the plane (77^, r]j), with 9 uniform over [0, 2ir), while if 
i = j we do nothing. 

Note that —C is a non-negative, bounded self-adjoint operator on L 2 (v). Any con- 
stant is an eigenfunction with eigenvalue and the spectral gap A = \(N) is defined 
as 

A W = inf *!££&, (1.2) 

/eL»: V{f 2 ) 

K/)=o 
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where v(f) stands for the expectation J fdv. M. Kac conjectured that A(A) stays 
bounded away from as N — > oo. This conjectured was first proved by Janvresse [7J, 
where a powerful recursive approach due to H.T. Yau was used. After that, in the 
beautiful paper [3], Carlen, Carvalho and Loss introduced a new recursive approach 
which allows to actually compute the value of A(A) for every N: 

AW~, N>2. (1.3) 

Around the same time, Maslen [9] derived formulae for all eigenvalues of C by means of 
harmonic analysis techniques. We refer to [3] for further background, motivation and 
references on Kac's conjecture. Our result below shows that the proof of (|1.3p can be 
somewhat simplified. In particular, we do not need any recursive analysis: in one step 
we go from X(N) to A(3) and the conclusion follows by direct computations in the case 
N = 3. 

Theorem 1.1. For any N ^ 3: 

A(A0 = (3A(3)-l)(l -£)+!. (1.4) 

In particular, U.3\) follows from Ji.^| j with A(3) = 5/12. 

The proof uses a well known equivalent characterization of the spectral gap as the 
largest constant A such that the inequality 

v{{Cff) > A i/ (/(-£)/), (1.5) 

holds for all / £ L 2 {v). The equivalence of (jl.2p and (jl.5p follows from elementary 
spectral theory. A similar approach has been exploited recently in [1] to obtain spectral 
gap bounds for a class of interacting particle systems. A proof of Theorem 11.11 is given 
at the end of the introduction. In the following sections we shall show that variants 
of the same method can be used to obtain spectral gap estimates for several models 
sharing some of the features of the Kac walk. The argument turns out to be especially 
powerful in non-homogeneous cases where other known methods are harder to apply 
because of the lack of permutation symmetry. In particular, for exclusion processes 
with site disorder we obtain a remarkable simplification of a spectral gap estimate 
proved by the author in [2] . The latter estimate is at the heart of recent results on the 
hydrodynamic limit of disordered lattice gases (6j [UJ . In the last section of this work 
we prove a new result which extends the spectral gap estimate to the case of colored 
particles. 

On the other hand, we point out that for some of the physically more relevant gen- 
eralizations of the Kac walk treated in [3] our argument will not necessarily yield sharp 
results. This is the case, for example, of the Kac model with momentum and energy 
conserving collisions discussed in the next section. The problem of the determination of 
the spectral gap for the latter model has been recently solved in [5], where the authors 
develop an interesting extension of the recursive scheme introduced in [3] . 

1.1. Proof of Theroem ll.il We start with some preliminaries. We write 

E b f( V ) = ±-J n f(R%ri)de, b = {i,j}, (1.6) 



for the binary average operator appearing in the definition of C. Note that E b is a 
projection which coincides with the ^-conditional expectation given the <r-algebra J- b 
generated by variables {rj£, £ £ b}. Thus, we rewrite the Markov generator as follows: 

£/fa) = ^£A/fa), (1-7) 

6 

where the sum runs over all ( „ ) unordered pairs b and for each such pair 

D b = E b - Id, E b f(ri) = u(f\F b ). (1.8) 

For each b, D b is a bounded self-adjoint operator in L 2 (u) satisfying D 2 = —D b . In 
particular, 

Hf(-C)g) = ^^2u[D b fD b g]. (1.9) 

6 

On the other hand we have 

«'((£/) 2 ) = ]^£ , WiV/]. (i.io) 

6,6' 

We are going to use the expressions (|1.9p and (ll.lOj) in (|1.5p to compute A(iV). We 
start with the lower bound. We write b ~ b' when two unordered pairs have at least one 
common vertex (including the case b = b'). Otherwise we write b 96 V '. We observe that 
if b 9^ b' , then E b and £/{,/ commute. Therefore, using D 2 = —D b and self-adjointness 

y[DbfD v f] = -v[{p v D b f){Dvf)] = v[(D b ,D b f) 2 } > , 6^6'. (1.11) 

It follows that 

^(( £ /) 2 ) ^ E KA./A,'/]. (1-12) 

6,6': 6-6' 

Unordered triples {z, j, A;} of distinct vertices are denoted by T (triangles). We say that 
b £ T if b = {i, j} and i, j £ T. Clearly, if b ~ b' and b ^ b' there is only one triangle T 
such that b, b' 6 T. We may therefore write 

E v[D b fD h ,f}= Y, v[D b fD b ,f\ + Y,v[{D b f) 2 ] 

6,6' : 6-6' 6,6' : 6~6' , 6^6' 6 

= E E ^DbfDyf] - E E ^(a,/) 2 ] + E ^(a,/) 2 ] • (Lis) 

T 6,6'GT T 6GT 6 

Since for every b there are exactly N — 2 triangles T such that b € T we see that 
E v[D b fD b ,f] 

6,6' : 6-6' 

= E E *WA,'/] - (iV - 3) E v[{D b ff] ■ (1-14) 

T 6,6'GT 6 

Let us now apply the inequality (II. 5p to a fixed triangle T. Let Tt denote the a- 
algebra generated by {rji, £ ^ T}. The conditional probability u[- \ Tt\ coincides with 
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the uniform probability measure on the sphere S 2 (t) in M 3 with radius 



Moreover, as noted in [1], it is not hard to show that the spectral gap of the Kac model 
does not depend on the radius of the sphere on which the walk is performed. Using 
(ll.5p on each triangle T, we therefore have 

l - Y v[D b fD v f | T T ] > A(3) Y v\W? | Ft] , (1.15) 

6,6'eT beT 

uniformly in rj £ S . Taking ^-expectation we can remove the conditioning on J~t 
in (|1.15p . Using this in (I1.14p gives 

Y v[D b fD v f] 

b,b> : b~b> 

> 3 A(3) (N-2)Y v[{D b ff] ~(N-3)Y *4(AJ) 2 ] 
b b 

= ((3A(3)-l)(iV-2) + l)^z,[( J D fe /) 2 ]. 

b 

From (jl.5p we conclude that X(N) is larger or equal than the right hand side of (|1.4p . 
It remains to show that this bound is attained for a given /. To this end, take 

N 

f N (v) = Y^ +const - ( L16 ) 

Let us first check that vlD^f^Dv In] = whenever b ^ b' , so that (|1.12p is an equality 
for f = f N . To this end, note that u[DbfNDb'fN} = ^[fNDbDb'fN] = 0. Indeed, if 
b = {i,j} and b' <ft b, then DyfNiv) depends on r]i, r]j only through rjf+rjj = l — Y^k^bVl 
and therefore DbDyfN = 0. 

Next, we need that (J1.15H is an equality as well. This requires checking that for 
N = 3, for any value of the conservation law ^«=i Vi = t > 0; /3 ls i U P to additive 
constants (that may depend on t), an eigenfunction of — C with eigenvalue A(3). For the 
solution of this 3-dimensional problem, as well as for the calculation of A(3) = 5/12, we 
refer to [H Section 3]. Once the estimates (|1.12|) and (|1.15|) can be turned into identities 
we see that all our bounds are saturated for the function (|1.16|) . This completes the 
proof. □ 

2. Homogeneous models 

The general setting can be described as follows. We consider a product space Q = 
X , where X, the single component space is a measurable space equipped with a 
probability measure \i. On Q we consider the product measure fj, N . Elements of Q will 
be denoted by r\ = (rji, ... , r/jy). Next, we take a measurable function £ : X —>■ R , for 
a given d ^ 1 , and we define the probability measure v = t>N,u on Q, as /i N conditioned 



on the event 

Sl NtU , := I n G fi : ££fa)=wl , (2-1) 

where w £ I is a given parameter. We interpret the constraint on CIn,u as a conser- 
vation law. 

In all the examples considered below there are no difficulties in defining the condi- 
tional probability v, therefore we do not attempt here at a justification of this setting 
in full generality but rather refer to the examples for full rigor. The crucial property 
of v that will be repeatedly used below is that, for any set of indices A, conditioned on 
the a-algebra Ta generated by variables rji, i £ A, v becomes the //-product law over 
7]j, j £ A, conditioned on the event 

We shall call this the non-interference property. In analogy with (II. 7|) we consider the 
binary collision Markov process described by the infinitesimal generator 

C f = jf E K/ I ft) - /] . (2-2) 

b 

where the sum runs over all ( 2 ) unordered pairs b = {i,j} and u[f | J~b\ is the u— 
conditional expectation of / given the variables rji, £ $. b. This defines a bounded 
self-adjoint operator on L 2 (u). Setting, as before, D^ = r/[-\ Fb] — Id, we see that 
D\ = — Db and the operator £ satisfies (jl.9j) . 

In principle, our arguments could be extended to more general processes. For in- 
stance, one can consider binary "collisions" which are not given by simple averages as 
in (|2.2p but by other mechanisms which still preserve reversibility (see e.g. the non uni- 
form models considered in [4]), or one could look at "collisions" which can involve more 
than two components at a time. However, we shall not investigate such extensions. 

Returning to our model (12. 2p . the spectral gap is defined as in (II. 2|) . Note that, by 
definition, A(2) = | always, since for N = 2 there is only one pair b and u((Dbf) 2 ) = 
v(f 2 ) for any / such that v{f) = 0. Here we shall assume that A(3) is independent of 
the choice of to. This is a strong assumption which holds only for special choices of the 
model. However, it does hold in the examples considered below. More general models 
are treated in the next section. 

Theorem 2.1. Suppose A (3) is independent ofu. Then, for any N Js 2: 

A(iV)M3A(3)-l)(l -£)+!. (2-3) 

If, in addition, there exists (p : X — ► R such that the function 

3 

Hvi^m^m) = ^Z^{m) , (2-4) 

2=1 

satisfies, for N = 3, Lfj, = — A(3)/3 + const., regardless of the value of Yli=i^,(Vi) 
(although the constant may depend on this value), then \2. 3\) can be turned into an 
identity for each N ^ 2. 
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Proof. We repeat the steps of the proof of Theorem 1 1. 11 We start from (jl.lOp and arrive 
at (|1.12p with the same commutation property used in (jl.lip . Indeed, this is a simple 
consequence of the non-interference property. The latter property also implies that the 
conditional probability v(- \ Ft) is nothing but v%^' with u/ = uj — YliiT^iVi)- Since 
A (3) is independent of the value u in the conservation law we may repeat the argument 
leading to (j!.15|) . This proves the lower bound (j2.3|) . As for the reverse direction, again 
the arguments given in the proof of Theorem 11.11 can be repeated line by line. □ 

Next, we examine some examples to which the theorem applies. 

2.1. Kac model. The model discussed in the introduction can be seen as a special 
case of our general setting, so that Theorem 11.11 becomes a special case of Theorem 
12.11 Here X = R, /x is the centered Gaussian measure with variance a 2 > and we 
take £(?7i) = ?yf (with d = 1). Then, for every uj > 0, vn,u is the uniform probability 
measure on the sphere of radius y/uj. Clearly, the choice of a 2 > is uninfluential in 
the determination of i/jv w . As we have seen in the introduction, this model satisfies the 
two main assumptions in Theorem 12.11 



2.2. "Flat" Kac model. Here X = M + , \i is the exponential law with parameter 
7 > 0. We take £,(Vi) = Vi (with d = 1). Then, for every ui > 0, independently of the 
choice of 7 > 0, vn,u is the uniform probability measure on the simplex Q,N tU . The 
binary collision process (|2.2|) associated to this setting does not appear explicitly in 
the literature, so we shall give more details on the computation of A(iV) in this case. 
Let us first check that X(N) is independent of uj, for any N. This is a consequence 
of the fact that C commutes with the unitary change of scale from Qjv,u to f2 jv,w' > for 
any ui, u/ > 0. Indeed, unui' is the image of vnu> under the map T : rj — > uj'tj/uj and if 
Mv) = f(Tv), then 

VN,u(fr I Fb)(ri) = VN,u'{f I Fb){Tri) , (2.5) 

for all r\ £ Qn,lu and for all pairs b. 

Next, we shall prove that A(3) = I and that the eigenfunction for iV = 3 is given by 

3 

h{m,m,m) = J2 r] i + const - ( 2 - 6 ) 

for any value of the constraint ^?.=i Vi = w (°f course, the constant will be given by 
— 3v3 >ul (r) 2 ), since we must have vs, w {f) = 0). From these facts and Theorem 12.11 we 
therefore obtain, for all N ^ 2: 

N + 1 

To solve the 3-dimensional problem, we observe that when iV = 3, then C + 1 coincides 
with the average operator P introduced in [4J. Therefore we can apply the general 
analysis of Section 2 in [JJ or equivalently that of Theorem 4.1 in [3]. The outcome is 
that 

A(3) > - min{2 + m , 2 - 2/i 2 } , (2.8) 



where the parameters fi\,fi2 are given by 

/ji = inf v(ip(rji)p(V2)) , ^2 = sup 7/(^(771)^(772)) (2.9) 

with tp chosen among all functions tp : X — ► R satisfying v(tp(r]i) 2 ) = 1 and v(tp(rji)) = 
0. Here v stands for 1/3^, but we have removed the subscripts for simplicity. As in 
(|2.5p one checks that the parameters fJ,i,/J>2 do not depend on u>. Write JCtp(a) = 
v[<p{t]2) \f]i = oc], oc ^ 0. This defines a self-adjoint Markov operator on L 2 {v\), where 
v\ is the marginal on 771 of v. In particular, the spectrum Sp(/C) of /C contains 1 (with 
eigenspace given by the constants). Then /ii,//2 are, respectively, the smallest and the 
largest value in Sp(/C)\{l}, as we see by writing v(ip{r\\)<p{r\2)) = v\tp{r\\)K,tp{ri\)\. This 
is now a one-dimensional problem and \i\,\ii can be computed as follows. To fix ideas 
we use the value u = 1 for the conservation law 771 + 772 + 773. In this case fi is the law 
on [0, 1] with density 2(1 — 771). Moreover, 

I rl-m 

£<p(m) = 1 / <p{m)dm , m e [o, i) • 

(K.cp(l) = </?(0)). In particular, yi(a) = a — 3 is an eigenfunction of /C with eigenvalue 
—1/2. Moreover, /C preserves the degree of polynomials so that if Q n denotes the space 
of all polynomials of degree d ^ n we have JCQ n C Q n . By induction we see that for 
each n ^ 1 the polynomial a n + q n -i(a), for a suitable g n _i G Q n -ii is an eigenfunction 
with eigenvalue /j n = ^~, <. , and it is orthogonal to Q n -i in L 2 {v\). Since the union 
of Q n , n ^ 1, is dense in l?{y\) this shows that there is a complete orthonormal set 
of eigenfunctions tp n , where tp n is a polynomial of degree n with eigenvalue /j n and 
Sp(/C) = {fi n , n = 0,1,...}. Therefore we can take /xi = — \ and \ii = | in the 
formula ([2T5]> . We conclude that A(3) ^ |. 

To end the proof we take / = rj{ + 7/2 + r/| and, using ^[77^ 1 772] = 5(772 — 2772 + 1) we 
compute 

4 

^/( ? ?) = "~g/W + const - 

Thus, A(3) = I and the eigenfunction is given by (|2.6p . Clearly, the unitary change of 
scale T introduced above does not alter the form of the eigenfunction so that all the 
hypothesis of Theorem 12. II apply and (12. 7\\ follows. 

2.3. Momentum and energy conserving collision model. . Here X = R 3 , and 

it is a centered 3-dimensional Gaussian law N(0, C) with covariance matrix C given 
by a multiple of the identity. Each coordinate r\i is a 3-dimensional velocity vector 
77", a = 1,2,3. We have d = 4 conservation laws, with £, a (r]i) = rjf, a = 1,2,3 
(momentum conservation) and £ (77^) = \r)i\ 2 = ^a=i( r ?f ) 2 (energy conservation). For 
any uj = (u) 1 , . . . , uj a ) € R 4 with w 4 > 0, v^,u is the uniform probability measure on the 
manifold £In,l> (whenever £In >uj / 0) . We refer to [3] for an explicit description of the 
probability measure and its main properties. It is still the case that A(3) is independent 
of the conservation law, see [3]. In particular, the lower bound (|2.3[) holds in this case. 
However, a computation of A(3) for this model shows that A(3) = | (see (12. lip below) 
and therefore the estimate becomes X(N) ^ jt which is rather poor. Indeed, it was 
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shown in [3] that the the spectral gap is bounded away from zero uniformly in N: 

inf X(N) >0. (2.10) 

Recently, by a very deep analysis of the Jacobi polynomials naturally associated to this 
model Carlen, Jeronimo and Loss [5] succeeded in computing X(N) exactly for every 

N: 

A(iV) = ^, N>3. (2.11) 

This shows our approach is too rough here. As we know from Theorem l2.1l the loss must 
come from the lack of the second property required in that theorem. It was shown in [5] 
that the eigenspace of X(N), for the choice uj = (0, 0, 0, 1), is spanned by the functions 

JV 



fN,a(v) = ^2\Vi\ 2 V? , a = 1,2,3. 



One cannot expect that the change of scale from ui to a/ transforms a linear combina- 
tion of /iv.a's into itself (up to multiplicative and additive constants), and the second 
property in Theorem 12.11 must fail here. 

Let us show that our approach can nevertheless be used to prove the weaker result 
(|2.10p without any recursive analysis. Namely, we prove that if A(4) > 1/4 then (|2.10p 
holds. The choice of triangles in the proof of Theorem 11.11 and Theorem 12.11 can be 
replaced by the choice of larger cliques (i.e. complete subgraphs) of the original complete 
graph. Namely, if in (11.130 we sum over cliques with 4 vertices instead of triangles we 
shall obtain the bound, for TV ^ 4: 

X(N) > (4A(4) - 1) Q - 1) + 1 , (2.12) 

instead of (12. 3|) . To see this, set for simplicity a^v '■= ^[DbfDy /], for a given / G L 2 (v), 
and recall that N 2 v[(Cf) 2 ] = Y^b V a b,b' ■ Denote by Q the cliques of 4 vertices and note 
that: for every b there are |(iV — 2)(N — 3) Q's such that b G Q; for any b ^ b' with 
b ~ b' there are (N — 3) Q's such that b, b' S Q; for any b, b' with b 96 b' there is only 
one Q such that b, b' £ Q. Then 

Y a w = j^^ Y Y ab > b ' 

b,b':b^b',b~b' Q b,b'£Q :b^b',b~b' 

= w^ Y \ Y ab > b> - Y ab > b> - 

Q {b,b'£Q b,b'eQ:b^b' 

^j^Y^-^Y^-j^Y Y *# 

Q b£Q Q b,b'eQ : b^b' 

= i(AT-2)(4A(4)- 1)^0^-^^ J2 a W- 

b 6,6' : b^b' 




Therefore 



/ 4 «&,&' = 2.^ «6,6' + 2^ ab ' b ' "*" Z^i ab ' b 
b,V b^b' : 6~6' b,b> : bjW b 

>i[(]V-2)(4A(4)-l) + 2]J>, 6 +(l 



^ ab,v 



6,6' : 6/6' 

>ipr-2)(4A(4)-l) + 2]J> 6 , 6) 

6 

where, in the last line, we have used (jl.lip . Since —Nu(fjCf) = Y2 b a b,b, this proves the 
claim (|2.12p . Note that this argument applies in the more general setting of Theorem 
12.11 and similar computations can in principle be carried out for any choice of cliques 
of k < N vertices. Therefore, to prove (|2.10p it suffices to prove A(4) > |. With the 
value A(4) = ± from (l2~TT1) this gives X(N) > ± + ^, for all N ^ 4. 

3. Non-homogeneous models 

We generalize the setting introduced above as follows. As before we take Q = X N 
but now each copy of X will be equipped with possibly distinct probability measures 
fii, i = 1,. . . ,N. Again we consider conservation laws as in (|2.ip . associated to a given 
function £ on X and a given parameter uj, and the probability v given by the product 
Hi X • • • X /xjy conditioned on rj £ £In,u- Note that the non-interference property still 
holds for this setting. The binary collision process is defined as in (12. 2D . Again, —C is 
a non-negative self-adjoint operator on L?{y) and we may define its spectral gap just 
as in (|1.2j) . This time, however, to keep track of the conservation law we shall write 
X(N,uj) instead of just X(N). Note that X(2,uj) = ^ always. We define 

X(N)=w£X(N,uj), (3.1) 

where the infimum ranges over the set of admissible values of the parameter uj. This 
set depends on the choice of the model and, as usual, we refer to the examples for fully 
rigorous formulations of the results. As a convention we may set X(N,uj) = +oo if 
uj is such that the measure v becomes a Dirac delta. Thanks to the non-interference 
property of v there is no difficulty in repeating the previous arguments to prove the 
following estimate. 

Theorem 3.1. For any N ^ 2 and any uj: 

A(A0M3A(3)-l)(l -£)+!. (3-2) 

Let us investigate some specific models. 

3.1. Non homogeneous Kac models. Consider the following non-homogeneous ver- 
sion of the "flat" Kac model introduced in Section [3 Take X = M + and /ij the proba- 
bility on M + with density 

— exp(-ry, +6j(r?j)), 

where 6j are bounded measurable functions and zi = J dxexp (— x + bi(x)) is the 
normalizing constant. We set B := supj |&i|oo- For this model, any value uj > of the 
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conservation law is allowed in the definition (|3.ip of \(N). We claim that for every 



e < | there is 5 > such that 



liminfA(iV) >e, (3.3) 

iV— >oo 



provided B < 5. Thanks to Theorem 13.11 to prove (|3.3f) it suffices to show that 
A(3, w) ^ |(1 — e(-B)) with some e(-B) — > as -B — ► 0, uniformly in w > 0, the value 
of the conservation law X^j=i r \i = u. To this end we shall use a standard perturbation 
argument. Let v denote the measure [i\ x /i 2 X ^(- 1 X)i=i % = w ) aim can ^o the same 
measure in the case b\ = b<i = 63 = 0. Thus, for any bounded measurable function g 
we have 

v{g) = ir- / <% / dr] 2 g{m,V2,u-Vi-V2)u(m,m), (3-4) 

^w JO JO 

where u(x,y) = e~ hl ^ c >~ h ^ y >~ bz ^ <j} ~ x ~ y > and Z w = f^ dx J^~ x dyu(x,y). Using 

e- 3B ^u( m , m )^e 3B , 
it is easily seen that, for any bounded / we have the bound between variances 

Var,(/)<e 6B Var V0 (/). (3.5) 

(Use g = (/ — i'o(f)) 2 m (|3.4p and the fact that z^(g) ^ Var;,(/)). The same reasoning 
shows that 

v(f(-C)f)>e- 10B Mf(-£o)f), (3.6) 

where Co is the generator corresponding to the choice b{ = 0. Indeed, 

6 t=i 



1 3 

-^^(Vai^/lrii)) , 



«=i 



(where Var t ,(/|r/j) denotes the variance of / w.r.t. v(-\rn)). For each j = 1,2,3 we 
have (as above) 



Var,(/ I m) > e- 4B Var V0 (/ | Vi ) , (3.7) 

uniformly in rj. A further comparison gives v (V&r u (f \i]i)) ^ e~ 10B uo (Var t , (/ 1 77^)) 
which implies (J3.6H . 

Recall that the spectral gap for £0 is equal to 4/9 regardless of the value of ui > 0. 
The previous estimates therefore imply that 

A(3,u,)^e- 16B . 

This proves our claim with e{B) = 1 — e _16B , from which (|3.3p follows. 

The same argument can be used to produce uniform lower bounds on the gap of 
non-homogeneous versions of the Kac walk on S and of the momentum and energy 
conserving collision model, under the assumption of small perturbations. In the latter 
model we need the argument in (I2.12p to obtain a uniform estimate inf/v A(iV) > 0. 
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3.2. Non uniform random permutations. We take X = {1, . . . , N} and, for each 

i = 1, . . . , N we consider a probability //j on X given by /^(r^ = j) = s —^, — , where 

bi : X — ► R are bounded functions, and Zj = X)i=i e~ bi ■ To model permutations we 
use A" conservation laws that will force all components of n to have distinct values: set 
d = N and define £ = (£ J )^ =1 with £ J (r/j) = l^^-j. Fixing cj = (1, 1, . . . , 1) we see 
that the set fijv.oj coincides with the set of AM permutations of A" letters. We define 
the probability v as usual by /ii x • • • /xjv(- | r/ G fijv,w)> Note that if the bias functions 
bi are all then v is simply the uniform probability measure over permutations. 

The binary collision is now a random transposition process. Note that only the value 
oj = (1, 1, . . . , 1) is considered for the conservation law so that X(N) = X(N,oj) for this 
model (there is no real infimum in (13. ip here). In the uniform case (bi = 0) a simple 
variant of the argument of Theorem 12.11 proves that X(N) = \ for every N ^ 2. 

While relaxation to equilibrium for the uniform case is well known, the non-uniform 
case is certainly less understood. Here we can show that if B = supj |6j|oo is sufficiently 
small we have inf jv X(N) > 0. More precisely, for every e < \ there exists 5 > such 
that B < 5 implies 

liminf X(N) ^e, 

N— >oo 

for all A^ ^ 2. To prove this we use exactly the same argument we have used to prove 
(|3.3|) . In particular, it suffices to show that A(3) — ► | as B — ► 0. The same perturbation 
argument will yield the desired estimate. To avoid repetitions we leave the details to 
the reader. 

3.3. Exclusion processes with site disorder. Here we consider a non-homogeneous 
version of what is sometimes called the Bernoulli-Laplace process. The inhomogeneous 
distribution models site impurities or site disorder. We take X = {0, 1}, /ij is a Bernoulli 
law with parameter pi & (0, 1), i.e. Hi(r]i = 1) = pi and [ii(j]i = 0) = 1 — pj. The value 
of r\i is interpreted as the presence (rji = 1) or absence (rn = 0) of a particle at the 
vertex i. The function £ is given by £(ry,) = r^ so that for any integer u 6 {0, 1, ... , N}, 
the set £In,lu denotes the configurations of uj particles over A^ vertices. The binary 
collision process (|2.2|) becomes nothing but the well known exclusion process on the 
complete graph {1, . . . , N}. This can be seen as follows. Given a pair b = {i,j} and a 
configuration r\ G &n,u, write u»y = uj — Yle&b %• Clearly u>ij G {0, 1, 2}. Observe that, 
if u>ij = 1 then u(f \ -Ff,)^) is given by 

« (1 -«> : /(„;l,0)+ „ "(' ,-"■„) T /(„;0,l) 



Pi(l-Pj)+Pj(l-Pi) ' Pi(l-Pj)+Pj0--Pi)' 

where, for simplicity we write explicitly the i-th and j'-th entries in f(rf) = f(n;rji,rjj). 
On the other hand in the case bJij G {0, 2} we have u(f \ ^b)(r/) = /(??). Setting 

cM = P^-P^-m) + Mi-PiMi-vj) (38) 

Pi{i -pj) +Pj{i -Pi) piii-pj) +Pj{i -Pi) 

we therefore obtain, for any 77 G fijv,w> 

"(/ 1 n){n) - m = cbiv) (f(v b ) - m) , (3.9) 
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where n b denotes the configuration in which rji and m have been exchanged. From (|3.9p 
wee see that C has the familiar form of the exclusion process. 

If we proceed by perturbative arguments (as in the previous two subsections) we 
would be able to prove a result of the form: if pi are (uniformly) sufficiently close to ^ 
then we have a uniform bound from below on the spectral gap. However, we shall prove 
here a much stronger result. We assume there exists e > such that the parameters pi 
satisfy 

e^pi^l-e, i = l,...,N. (3.10) 

The minimal spectral gap \(N) is defined as usual by (|3.ip . where the infimum ranges 
over all uj 6 {0, 1, ... , N}, with the convention that X(N, 0) = \(N, N) = oo. Under the 
same assumption the following theorem was proved in [2 J by means of rather technical 
local limit theorem estimates. It is surprising that the simple argument of Theorem 13. II 
allows a straightforward proof. The uniform spectral gap bound below is an important 
step in the recent works [El II] establishing hydrodynamic limits for exclusion processes 
with disorder. 

Theorem 3.2. Assume 113. 10\) for some e > 0. Then, there exists c e > such that for 
all N^2 

\{N)>c e . (3.11) 

Proof. Thanks to Theorem 13. H all we have to do is prove that 

A(3) >- + <%, (3.12) 

for some c £ > 0. We fix three vertices i = 1,2,3, with their occupation probabilities 
Pi satisfying (|3.1U|) and with J2 i=1 rji = uj. We may assume that uj = 1, i.e. there is 
one particle. Indeed if there are two we may look at occupied vertices as empty and 
vice- versa, if there is none (or three) the measure is a Dirac delta and by our convention 
(see discussion after (|3.ip ) the estimate A(3, w) ^ | + c £ becomes obvious. Since there 
is one particle we shall call x, respectively y, the probability that the particle is at 
i = 1, respectively i = 2. We set z = 1 — x — y for the probability that the particle is 
at i = 3. Note that, thanks to (|3.10p . x,y,z are all bounded away from and 1. For 
instance, 

Pl(l -P2)(l -Pi) 

x = 

Pl(l -P2)(l -P3) + (1 -Pl)P2(l -P3) + (1 -Pl)(l ~P2)P3 ' 

It is easily checked that the process generated by C on our three sites becomes a 3-state 
Markov chain with the 3x3 transition matrix P = C + Id given by 

/l + ^ + ^ L - - JL - -2- \ 

i / x+y x+z x+y x+z \ 

r ~ o I x+y ' x+y ~ y+z y+z J ■ {O.iO) 

\ -Z- -U- 1 + ^ + ^/ 

\ x+z y+z x+z y+z / 

We need to estimate the eigenvalues of P. Clearly, one eigenvalue is 1. From (13.131) we 

see that Tr(P) = 2. Therefore the other two eigenvalues must satisfy Ai = Tr(P) — 1 — 

A 2 = 1 - A 2 . Note that 

A(3, w) = min{l — Ai, 1 — A 2 } • 

To estimate Xi,i = 1,2 we compute the determinant of P. The next lemma shows 

that det(P) > |. Therefore, for both i = 1,2 we have Aj(l — Aj) > |. This implies 
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|Aj — 2 1 < g- I n particular, it shows that Aj < |. In conclusion: A (3) > „- as claimed. 
(Note that Aj = ^ would be the value in the homogeneous case pi = g-) 

Lemma 3.3. 

det ' p »-K i+ a-^-;)(i-J <3j4) 

Proof. Set T = 3P. Also, set 

x x y 

* = — — > /?= — — > 7=— — • 

x + z x + y y + -2 

From (J3.13P we compute 

det(r) = (l + l3 + 5)(6-3f3-25 + f35-~f5 + j(3) 

+ (P - I) (3(3 - j3 6 - 7 P - 6 + i 6) 

+ (l-6)([3 6-i 6 + i (3-26) . 
Simplifying we arrive at 

det(r) = 6 + 35-3^5-37(5 + 3/37. 
Rewriting this in terms of the probabilities x, y, z we obtain 

det(r) = 6 + - J, Xy \ n v- ^ 3 - 15 ) 

(1 - x) (1 - y) (1 - z) 

This implies (l3~T31) . D 

3.4. Colored exclusion processes with site disorder. A natural generalization of 
the previous model is a system where particles can be of several different kinds - or 
colors. Namely, suppose there are m colors and let each particle be painted with one 
of the available colors. Configurations of colored particles are denoted by r] £ Q := 
{0,1,..., m} with the interpretation that r\i = means that i is empty, while r}% = k, 
k £ {1, . . . , m} means that i is occupied by a particle with color k. Thus, the single 
state space X is {0, ...,m}. The conservation laws are given by £ fc (f?j) = l{^=k}, 
k = I, ... ,m and the vector u> = (u>i, . . . , u> m ), where u>i are non-negative integers such 
that YlT=i ^k ^ N. Thus, the set Qn >u1 denotes the configurations of particles over A 
vertices, with lo^ particles with color k. We shall use the notation ifti = lr^. ^ u so that 
the variables tp £ {0, 1}^ denote the configuration of occupied sites. Let pi , i = 1, . . . A 
be given parameters satisfying (|3.10p . Let /ij denote the probability on X such that 

Viim = k ) = V" (Pi l {k > 1} + U - Pi) l {k=o}) , 

where Zi = (m — l)pi + 1 is the normalization. In particular, w.r.t. /x, the occupation 
variable ipi is a Bernoulli random variable with parameter mpi/Zi. We call, as usual, 
v the product \i\ x •••[in conditioned on f2jv>- To avoid degenerate cases we take 
1 ^ ^!k=\ LO k ^ A — 1. We are interested in the following exclusion-type dynamics. 
For any configuration 7] and any edge b = {i,j} we write rj b for the configuration 
where the variables rji,r]j have been exchanged. For every 7 ^ we define the Markov 
generator by 

W fa) = J3 E C ^V) (/fa 6 ) ~ I(V)) , (3.16) 

b 
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where the rates are given, for b = (i,j), by 

cjiv) = cb{ip) + 2 l {i>i=i> 3 } ( 3 - 17 ) 

where Cb are the functions defined in (j3.8fl . but now evaluated at the occupation vari- 
ables tp = ip{rf) defined above. It is easily checked that the rates (|3.17p are reversible 
w.r.t. v. for any n and any b 

u( V )cJ( V ) = u(r, b )cJ( V b ). (3.18) 

The latter statement is equivalent to self-adjointness of £ 7 in l?{y). Moreover, 



,(/£ 7 /) = ^.[c^(/ b -/) 2 ] , (3.19) 



where / (77) := f(jj). If 1 ^ SfcLi w fc ^ N — 1 the Markov chain generated by £ 7 is 
irreducible. We shall consider the cases 7 = and 7 = 1. The difference between Cq 
and L\ is that in L\ we have added the possibility of "stirring" between particles, i.e. 
exchange of positions of particles of different colors. The case 7 = 1 coincides then with 
the usual binary collision dynamics given by local averages (|2.2p . On the other hand, 
the case 7 = is a true exclusion process, with particles jumping only to empty sites. 
The addition of stirring can result in a faster relaxation to the equilibrium distribution 
z/, or equivalently in a larger spectral gap. Note, however, that for the case m = 1 there 
is no difference: in this case £ 7 = Cq for all 7. The case m = 1, of course, is the one 
analyzed in Theorem 13.21 From now on we take m > 1. 

We denote by A 7 (iV, uj) the spectral gap of the generator £ 7 . The following theorem 
shows that in the case 7 = 1 we have a uniform lower bound \\{N,lo) ^ c £ as in 
Theorem 13.21 and, in the case 7 = we have Ao(A r , uj) ^ c £ (l — p), where p is the global 
density: 



in 



fc=i 

The slow-down in the limit p — * 1 is natural in view of the absence of stirring. Moreover, 
we shall show that, up to a constant the reverse inequality \q(N,uj) ^ C(l — p) holds 
as well for some choices of uj, see the remark after the end of the proof. Similar results 
had been obtained in [10J in the homogeneous case pi = 7^. 

Theorem 3.4. Assume \3. 1 0\) for some e > 0. For m > 1, 7 = 1, there exists c £ > 
such that for any N ^ 2 and any uj such that < p < 1: 

X 1 (N,uj)^c £ . (3.21) 

For m > 1, 7 = 0, there exists c e > such that for any N ^ 2 and any uj with 
< p < 1: 

X (N,uj)^c £ (l-p). (3.22) 

Proof. We start with some preliminary facts. Consider functions / of the occupation 
variables i/j = {ipi} only. Let TCq denote the space of all such functions and observe that 
C^Tio C TCo, i.e. Tio is invariant, for both 7 = 0, 1. This follows from the fact that the 
only dependence of the rates on the configuration n is through the variables {ipi}, see 
(|3.17p . In particular, under the generator £ 7 , the variables {ipi} evolve as the Markov 
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chain of the case m = 1. Therefore, the spectral gap estimate of Theorem 13.21 applies 
to functions in TCq, for both 7 = 0, 1. 

For any / G L 2 (u) we may write / = /o + /o", where /o = u(f \ T^) 6 "Ho and 
fo~ = f — fo£ "Ho- Here ^ denotes the <r-algebra generated by the functions ipi and 
H.q is the orthogonal complement of Ho, i.e. the space of / such that v{fg) = for all 
g G Ho- Since £ 7 leaves Ho invariant, by self-adjointness it follows that C^Hq C Hq. 
In conclusion 

Moreover, Var I/ (/) = Var^(/o) + Var t ,(/o L ). From these simple observations it is clear 
that the constant A 7 (iV, u) satisfies 

A 7 (iV,u/)^ miii{A°(.W,u;),A^-(i\r,w)}, (3.23) 

where X^(N,u),X^r(N,u>) are the constants obtained in the variational principle (|1.2|) 
- applied to £ 7 - by restricting to / G 7^o and / G Hq respectively. From Theorem 
13.21 we know that 

A 7 (iV, W )^ C£ , (3.24) 

for some c e , for both 7 = 0, 1. To prove Theorem 13.41 we thus have to estimate from 
below the constant Ar(JV, u). 

3.4.1. The case 7 = 1. If / G Hq, then we must have u(f \ T^) = 0. Therefore 

Var„(/) = i/[Var(/ | ^)] + Var„ (u(f \ T^)) = i/[Var(/ | ^)] . (3.25) 

Here Var(- 1 .F^,) denotes variance w.r.t. v{- \ J-f), the conditional probability given the 
color-blind configuration tp, and we have used the standard decomposition of variance 
under conditioning. Let S = S(r)) := {i G {1,...,N} : ipi = 1} denote the set 
of occupied sites. Clearly \S\ = pN , the total number of particles. Observe that 
once i/j is known then the distribution of the variables r\ is given by rji = for % ^ 
S (deterministically) and rji G {1, ...,m} on S uniformly with the constraint that 
J2ies ^-{vi=k} = °°ki k G {1, . . . , m}. In particular there is no inhomogeneity once the 
set S (or, equivalently the configuration ip) is given. Therefore Var(/ | T^) can be 
estimated for every r\ with the well known bound for random transpositions without 
disorder (see e.g. [2]): 

Var(/ I ^) < 41W I> t( V ^'/) 2 I ?*] • ( 3 - 26 ) 

Here we use the notation Vijf(rj) = f{i] b ) — /(??), b = {i,j}, for the exchange gradient. 
Averaging with v and using (|3.25p we obtain (since \S\ = pN, deterministically) 

Vai„(/) ^ j^v v ( Y, ( V ^/) 2 I • ( 3 - 27 ) 

9 \ijeS J 

Suppose first p ^ i. Then (|3.27p . (13.191) and a uniform lower bound on the rates q, 
imply 

V«,(/)<2l E^[( V ^) 2 ] <C^(/(-A)/), 
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for some constant C £ < oo, with the sum now extending to all pairs i,j. This shows 
that, with c £ = l/C £ : 



\i{N,uj)>c £ , p>\. 



(3.28) 



Note that here we have used 7 = 1 (if 7 = there is no uniform lower bound on the 
rates q,). 

We turn to the proof in the case p ^ i . We rewrite (|3.27|) as 



Va ^ (/) ^ AjN ?> [ (V ^ /)2 l ^ l ^~\ 



1 -j 






Ap(l-p)N 2 



(3.29) 



where we use the identity (1 — p)N = N — J2k=i w & = Yle ^-{££S} f° r the number of 
empty sites. Let r\ E £In,u be fixed. Pick i,j E S(rj) and write r/ lJ for the exchanged 
configuration. Observe that for any £ ^ S(rj) we can write 



rf' 3 = \W 



1 -j 



Therefore 



vvf(v) = mv id )-m] 



'■j 



Vjtfllv 1 '* )+Vijf[v ht : +v«/(t/). 



ii 



Each term in the sum appearing in (|3.29|) can then be estimated by 



< 3v 



VjtflW' 



':) 



A,t 



(3.30) 



v«/ k + (yufr } i {jes} i 0G5} i W s } 



Next, we claim that there exists C\ = C\(e) < 00 such that 



1 -j 






1 {«es}l{jeS}l{^S} 



(V^/(r/)) l{ ie5 }l{^ S }l^ e5 } 



(3.31) 



17 



Note the change of the indicator functions in (|3.3ip . To prove (|3.3ip we write, with 
the change of variables (p := rf' and := ip^'i 



l/,e<nli 






Lli 



{iGS} i Oes}i{^S} 

2 



;.c 



'■j 



1 {iGS(^)}l0'GS(^)}l{^S(^)} 



EK^)^/(^)) 2 %%)}1 



b'eS( v .)}l{<GS( v .)} 



E^((^ J ) 4 ^) (Vj£/(/3)) 2 l{ie5(/3)}l{i^S08)}l{^e5(/3)} 



^ [Xtj,^ (Vjif) l{igS}l{j^S} 1 {^eS} 



where 



XnWfa) := 



v({rfjy> 1 ) 



u(rj) 

is the change of measure. Since the variables p, defining our measure satisfy (I3.10P it 
is not hard to check that Xi,j,i(v) ^ C e uniformly for some constant C £ . This proves 

Moreover, in a similar way one proves that there is a constant C2 = ^(e) < 00 such 
that 

2 



(Vi;/ 



;.f 



'/' 



l.ric.<n.l.r..-c.<?il 



{iGS}i{j"GS}i{^S} 



^<? 2 ^ 



(3.32) 



(Vij/(r?)) l{^5}l{ ie5 }l{ te5 } 

From (13.291) and (j3.30p . using ()3.3ip and (13.321) we obtain for a suitable constant C3: 
C 3 



V ^ (/) ^ ,fl- 3 ^ E {^ t( V ^) 2 ^S}!^^}! 



{^eS}J 



p(l-p)*\,, 

+ v [(Vijf) 2 l{^S}l{iGS}l{«eS}] + " [(V^/) 2 l {ie5} l {j6 5 } l W 5 } ] J . (3.33) 

Since J2% ^{ieS} = P-^> setting C4 = 3 C3 we see that (|3.33|) can be written as 

C 4 



Var "^ < 7T37wE ,/ [( v «/) 2l {«65}l^ 



5} J 



(3.34) 



(1 - p)N 

h3 

Using p ^ ~ and (|3.19|) we see that 

Var„(/)<Cz/(/(-£i)/), 

for some constant C = C(e) < 00 and for all / £ Wq ■ Therefore we have proved 
that Xi(N,uj) > c e . Together with (pT23j) . (pTHj) and (pTZHj) we see that Xi(N,p) is 
uniformly bounded from below and the claim (I3.21J) holds. This proves Theorem 13.41 in 
the case 7=1. 
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3.5. The case 7 = 0. Here we cannot use the argument leading to (|3.28[) above. 
However, we can use the argument giving (|3.34p without modification (it never used 
the fact that 7=1). In particular, (|3.34p holds for any < p < 1 and any / G T~Cq ■ 
Now, observe that for any edge b = {i,j} such that i £ S and j ^ S the rate c° defined 
in (|3.17|) is uniformly bounded away from zero with constants only depending on the 
e from (|3.1U|) (this follows from the fact that for such cases either ipi(l — ipj) = 1 or 



±. 
• 



ipj(l — ipi) = 1). Therefore, there exists C = C(e) < 00 such that for any / 6 TL 

Var„(/) < — ^— i/(/(-A,)/) • (3.35) 

(I-P) 

This proves that \£(N,uj) > c(l - p). From (l3~23l) and (pOlj) we see that A (A r ,w) 
satisfies the claim (|3.22p . This proves Theorem 13.41 in the case 7 = 0. □ 

Let us briefly address the issue of comparable upper bounds on spectral gaps. For 
example, to prove that Xo(N, to) = 0(1 — p), as p — ► 1 one can consider the case of 
m = 2 colors, with u)\ = 002 such that u\ + lo<i = pN with p — > 1. Then choose / in 
the variational principle (|1.2|) as the indicator function of the event that a given vertex 
i is occupied by a particle of color 1. The variance of / is of order 1. On the other 
hand, since the number of empty sites that can be used to change the value of rji is 
(1 — p)N, it is not hard to show that v(f(—Co)f) is of order (1 — p). It follows that 
Ao(iV, to) = 0(1 — p). Of course, if 7 = 1 then v(f(—C\)f) is of order 1 and the gap 
does not vanish as p — ► 1 in that case. 

Finally, we point out an interesting problem concerning local versions of the exclusion 
dynamics described by (|3.19p . The local dynamics is obtained by summing over pairs 
b which are given by the edges of a small subgraph of the complete graph, such as e.g. 
the box of side L ~ JV 1 '^ in a d-dimensional grid: A^ = [1,-L] d n Z d . In the latter 
cases one expects the gap to scale as L~ 2 . In the case m = 1 there is a nice argument 
in [11| Lemma 5.1 and Lemma 5.2] which allows to derive such an estimate from the 
complete-graph bound (|3.11|) . On the other hand, the case m > 1 with 7 = seems to 
be much more complicated and we are not aware of any result in this direction, except 
for the homogeneous case considered in [IT 
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